Typeset with jpsj2.cls <ver.l.2> 



Full Paper 



Crystalline-Electric-Field Effect on the Resistivity 

of 

Ce-based Heavy Fermion Systems 

Yasutaka Nishida*, Atsushi TSURUTA and Kazumasa Miyake 

Department of Materials Engineering Science, 
Graduate School of Engineering Science, Osaka University, 
Toyonaka, Osaka 560-8531 
(Received February 6, 2008) 

The behavior of the resistivity of Ce-based heavy fermion systems is studied using a 1/ A^- 
expansion method a la Nagoya, where N is the spin-orbital degeneracy of f-electrons. The 
1/ A^'-expansion is performed in terms of the auxihary particles, and a strict requirement 
of the local constraints is fulfilled for each order of The physical quantities can be 

calculated over the entire temperature range by solving the coupled Dyson equations for 
the Green functions self-consistently at each temperature. This 1/A^-expansion method is 
known to provide asymptotically exact results for the behavior of physical quantities in 
both low- and high-energy regions when it is applied to a single orbital periodic Anderson 
model (PAM). On the basis of a generalized PAM including crystalline-electric-field splitting 
with a single conduction band, the pressure dependence of the resistivity p is calculated by 
parameterizing the effect of pressure as the variation of the hybridization parameter between 
the conduction electrons and f-electrons. The main result of the present study is that the 
double-peak structure of the T-dcpendence of p{T) is shown to merge into a single-peak 
structure with increasing pressure. 

KEYWORDS: Ce-based heavy fermions, electrical resistivity, 1/A^-expansion method, crys- 
talline electric field 



1. Introduction 

The resistivity of heavy fermions such as Ce-based compounds (Kondo-lattice systems) 
changes markedly with decreasing temperature. In the high-temperature region the conduction 
electrons are decoupled from the f-electrons, which behave as localized moments, exhibiting the 
Kondo effect that results in the logarithmic temperature dependence of the resistivity. On the 
other hand, below the characteristic temperature Eq (corresponding to the Kondo temperature 
in an impurity model) conduction electrons are scattered coherently by the f-electrons, and a 
Fermi-liquid state with quasi-particles of increased mass is formed. The detailed T dependence 
of the resistivity p{T) around T = Eq can be roughly classified into two cases. In the first case, 
p{T) has a single-peak structure as observed in CeCue-^^ In the second case, p{T) has a double- 
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peak structure arising from the crystalline-electric-field (CEF) effect, as observed in CeAl2,^^ 



CeCu2,'^^ CeAu2Si2^^ CeCu2Si2,^^ and CeCu2Ge2,^'* for example. Experimental studies of the 
effect of pressure on Ce-based compounds have revealed the fact that the double peak in the 
resistivity curve arising from the CEF effect tends to merge into a single peak as the pressure 
is increased. In the 1980's the effect of CEF on the impurity Kondo effect was investigated 



heavy fermions (lattice systems) has not been adequately investigated to date. In this work, 
we study the influence of the CEF effect on the T-dependence of p of Ce-based heavy fermions, 
and clarify the origin of the double-peak structure of p{T). 

To this end, we investigate an infinite-f7 generalized periodic Anderson model (PAM), in- 
cluding the CEF splitting of f-electrons, by using a 1/ A^-expansion method d la Nagoya, where 

is the spin-orbital degeneracy of correlated electrons.^' ^'^^ This 1/A^-expansion method is 
performed in terms of the auxiliary particles together with strict requirement of the local con- 
straints in each order of 1 /N so as to rigorously take into account the effects of the infiniteness 
of U . As a result, the theory is free from an artificial transition of slave-boson condensation 
and gives correct behaviors of physical quantities in both low- and high-energy regions. Apply- 
ing this method, we discuss the pressure effect of p{T), and show that the result is consistent 
with the tendency observed in Ce-based heavy fermion systems. 

The paper is organized as follows. In §2, we introduce the model Hamiltonian, a generalized 
PAM, and the concept of the 1/A^-expansion method d la Nagoya. In §3, we show the result at 
T = (the case of absolute zero temperature). In §4, we show results at a finite temperature by 
using a self-consistent 1/A^-expansion beyond the low-temperature approximation. Finally, we 
show the effect of pressure on the resistivity p{T) in §5. Three appendices are given detailing 
the calculation and algorithm of the 1/A^-expansion method. 

2. Model and Formal Preliminaries 

2.1 Infinite-U generalized periodic Anderson model 

We start with an infinite-C/ periodic Anderson model (PAM) with an f-electron in a man- 
ifold of J = 5/2,11) ^j^ipj^ 

is a simple and realistic model for the Ce-based heavy fermion 
systems: Our model Hamiltonian is given by 



theoretically by various methods.'^' ^'^'^'^^^ However, to our knowledge, the effect of CEF in 



H =Hc + Hf + Hcf 



(2.1) 




(2.2) 




(2.3) 




(2.4) 
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|±3> = b|=F5/2> + a|±3/2> 



a +b =1 



Fig. 1. CEF level scheme of Ce ion in tetragonal symmetry. 



where c|^^ is the creation operator for the conduction electron with wave vector k and spin 
(7, and ffkcr is the energy of the conduction electron with wave vector k and spin a. 
denotes the number of the lattice sites, and i stands for the site index. Here, we introduce the 
slave-boson, which represents the /'^-state at each lattice point, following Coleman:^^) h\ is 
the creation operator for the slave-boson at the i-th site. /Jp is the creation operator for the 
pseudo-fermion which represents the /^-state with the CEF level |r) = | ±1), | ±2), | ±3). Vkr^ 
stands for the mixing between the conduction electron with k, a and the f-electron with the 
CEF state of F. We discuss the property of the generalized PAM introduced above assuming 
under a tetragonal symmetry of the crystal field. Following Kontani et al.,^^^ the mixing Vkro- 
is expressed as 



M 



OtmV\^Mu, 



V]iMa = Vq 





3 1 -V ^-^ri^^y 



(2.5) 



(2.6) 



where M is the z-component of the angular momentum J of an f-electron: M = {J = 
5/2), and Yj^{Q.k) is the spherical harmonic function. In the tetragonal symmetry, we can 
parameterize each CEF-level as follows: (See Fig. 1.) 



(2.7) 



ej is measured from the chemical potential £/ = — and e^p is defined as the lowest 
atomic level of f-electrons measured from the center of the conduction electron band. The 
coefficient OrM in eq. (2.5) is a 6 x 6 orthogonal matrix, which is given as follows: 
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Fig. 2. Feynman diagrams: A solid line denotes the conduction electron propagator, a wavy line the 
slave-boson propagator, and a dashed line the pseudo-fermion propagator. 



Hereafter, all the energies are measured from the chemical potential ^ = 0. 

Next, we recapitulate the 1/A^-expansion method d la Nagoya. It is straightforward to 
extend the work by Ono et al. for a single-band PAM as detailed in ref. 10. In order to 
ensure equivalence between the present model Hamiltonian eqs. (2.1)~(2.4) and the original 
infinite-C/ generalized PAM, we must treat the problem under the local constraints as 

Q^ = h\hi + J2 flvhr = 1- (2.9) 
r 

The expectation value of an operator O under the local constraint eq. (2.9) is given as 

(0)= hm {6T\Q,:)x/{T\Qr)x, (2.10) 

where {A)\ is calculated in the grand canonical ensemble for the Hamiltonian H\: 

{A)x = TrJe-'^-^Mj/Tr^e-^^^] , (2.11) 

Hx = H + Y,>^iQ^■ (2.12) 

i 

The single-particle Green functions Gko- for the conduction electron, Bi for the slave boson 
and Fjr for the pseudo-fermion are essential ingredients. Their unperturbed forms are given 
as follows: 

Gi^iLOn) = (iiOn - ek^)-\ (2.13) 

B^{iUn) = {wn-Xi)-\ (2.14) 

F^^iojn) = (iiOn -Xi- Eir)-\ (2.15) 

where tOn = (2n+l)7rT, and Vn = 2mrT with n being an integer. Then, the Feynman diagrams 
are illustrated in Fig. 2, and those for the c/-mixing vertices are illustrated in Fig. 3. The 
effect of the interaction, eq. (2.4), can be formally incorporated into the Green function as 
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self-energies as follows: 

Gkai^^n) = (i^n " Eka " ^ka{^^^n))~^ , (2.16) 

Bi{iu^) = iivn -Xi- U,{wn))-\ (2.17) 

FiriiuJn) = {iujn - Eir - Xi - sfp(ic^„))-^ (2.18) 

The self-energies Sko-i Hj, and S^p are calculated by an extended 1/ A^-expansion d la Nagoya, 
as discussed below. 

2.2 1/N-expansion d la Nagoya for SU(N)-PAM 

We briefly review the past work performed by Ono et al. as detailed in refs. 9, 10. The 
conventional 1/A^-expansion method for SU(A^)-PAM, in which the conduction band is pre- 
pared for each component of pseudo-fermion and the hybridization depends only on k, is 
performed by dividing the k-space and spin-space of the conduction electron into N subspace 
and keeping the total degrees of freedom of the conduction electron to be the same as those 
of the non-interacting conduction band: 

{T, k] + {i, A;} = {fell + {^2} + • • • {km] + ■■■ + {feiv}, (2.19) 

m=J ^ 

where Nl stands for the number of total lattice sites, A'^^: = X^i ^- assume that the km 
state of the conduction electron hybridizes only with the m-th state of the f-electron. Then, 
the density of states (DOS) of the m-th subspace of the conduction electron is 1/A^ times as 
small as the total DOS, so that a summation over km gives a factor of \/N. Therefore, we 
can classify the Feynman diagrams giving the self-energies of the Green functions in terms of 
1/A^. This method may be justified on the basis of the following intuitive argument. 

The state of the conduction electron with momentum k and spin a is expanded around 
the origin, where an f-electron is located, into that of the f-electron with the state of angular 
momentum state \m): 

Cka = '^'^{cr,'k\m,km)ck^- (2.21) 

m km 

Namely, only a restricted part of the degrees of freedom k and a of the conduction electron 
near the Fermi level hybridizes with f-electron in the state |m). If the f-electron state is 
specified by N different quantum numbers, the degrees of freedom of conduction electrons 
which hybridize with the f-electron with the state \m) are equal to the roughly 1/A^ times of 
the total degrees of freedom of the conduction electron, so that the effective DOS of a relevant 
conduction electron is reduced by 1/A^ times. In the lattice system, hybridization between 
conduction electrons with k and a and the f-electron with |m) at the i-th site is given in 
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Fig. 3. Fcynman diagrams for cf- hybridization: A solid line denotes the conduction electron propa- 
gator, a wavy line the slave-boson propagator, and a dashed line the pseudo-fermion propagator. 
These propagators are connected by the cf-hybridization vertex due to Hcf, eq. (2.4). 



terms of the matrix element in (2.21) as 

{a,'k\km,m;i) = e~'^''^' {a,'k\m, km) ■ 
Therefore, the Hamiltonian Hcf for hybridization is given as follows: 
Hcf =^^(^(cr,k|i7|m;i)4^/im -hh.c. 

i,m cTjk 



^ ^ '^(^{rn',krn'\(T,k){a,k\H\rn-,i)cl^Jirn + h.c^, 

i,m m',k^i (T,k 

X] [{^^'^km'\H\m;i)cl^Jim + h..c^, 



i.m m'k^i 



i,m m'k^i 



i,m kn 

Similarly, the Hamiltonian He for the kinetic term is given as follows: 

o-.k cr'.k' 

= E E ^^i'^^'^rnW,k){a,k\H\a'y){a',k'\m',krn')cl^Ck^,, 

m,kmm' ,k^i (T,k g-'jk' 
m,km m',k^i 



2.22) 

2.23) 
2.24) 
2.25) 
2.26) 
2.27) 

2.28) 
2.29) 
2.30) 
2.31) 
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Here, we have replaced the matrix elements {m' ,km'\H\m;i) in eq. (2.25) with 
VSm,m' e"''""™'^' / a/ Nl , and we have replaced {m, km\H\m' , km') in eq. (2.30) with 
^km^m,m'^km,k^/ ■ Then, the rareness of the DOS of the conduction electron for each channel 
m is inherited in the lattice system. In the lattice system, strictly speaking, the conduction 
electrons described by Cfc,^'s mix with each other so that a classification in terms of the index 
m becomes ambiguous. So the approximation in eqs. (2.26) and (2.31) for the conduction elec- 
trons should be regarded as an approximation that neglects the hybridization of conduction 
electrons with different indices m. Nevertheless, this approximate treatment has been shown 
to provide reasonable results for numerous strongly correlated electron systems even for the 
case iV = 2.1^) 

However, an extension of the conventional SU(A^)-PAM enabling inclusion of CEF split- 
ting A causes a few problems: e.g., the pole of Bi{v), Eq which corresponds to the Kondo 
temperature Tk, results in contradiction, i.e., \ + £f — Eq > A-|-e/, if A is larger than almost 
twice Tk for A = 0, i.e., without CEF splitting. In a word, the energy level of the coherent 
state becomes higher than that of the bare f-electron, and the 1/ A^-expansion method fails. 
For example, in the case where four degenerate states of the f-electron, i.e., = 4, are di- 

(2) 

vided into two Kramers doublets by CEF splitting A, the correct Kondo temperature 
for = 2 is not obtained in the limit A ^ od. One of the origins of this difficulty may be 
that the hybridization between conduction electrons with different quantum numbers m is 
discarded, so that the two components of the conduction electrons are independent of each 
other. To resolve this problem, we reconsider a generalized PAM in which the conduction 
electron, specified by k and a, mixes with all the states of f-electrons specified by m. 

2.3 Extended 1/N-expansion method for generalized PAM 

We now apply the 1/A^-expansion method to the generalized PAM introduced in §2.1. 
Namely, we extend a conventional 1/A^-expansion method in such a way that the condition 
of eq. (2.20) is replaced by the following conditions. 

l-Y^l = l = Oil/N), (2.32) 
Y^l = 2 = 0{N), (2.33) 

a 

]^EEl = 2 = 0(l). (2.34) 

Then, the degeneracy A^ is given by just that of the real spin of the conduction electrons: 
i.e., N = 2. We use this condition hereafter in the present paper. As shown below, the rule 
of classification of the Feynman diagrams with respect to 1/A^ in the SU(A^)-PAM model can 
be generalized in the present model that has only spin degeneracy N = 2. The validity of 
using the 1/A^-expansion method d la Nagoya in the case A^ = 2 is supported by the results 
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(3 a 

c l'„ k' 




Fig. 4. Example of diagrams showing contributions to the self-energy: the lowest order of the self- 
energy for the slave boson (a) and the next order examples of the self-energy for the conduction 
electron (b,c,d) in power of (1/-/V). 



of previous studies that treated a single-band PAM with N = 2 and gave quahtatively correct 
physical behavior in both the high- and low-energy regions. In particular, it has been shown 
that higher order corrections in do not change the results qualitatively.^^) On the basis of 
this observation, we use this 1/A^-expansion as a method of taking into account the correlation 
effect of f-electrons in the present model. 

The self-energies of the conduction electron Ska-, in eq. (2.16), the slave-boson Ilj, in eq. 
(2.17), and the pseudo-fermion sjp, in eq. (2.18), are determined for each order of 1/A^. To 
classify the diagrams by each order of (1/A^), we use the condition in eqs. (2.32) (2.34). 
We illustrate several examples of the self energy in Fig. 4. In our extended formulation of the 
1/A^-expansion, the diagram of Fig. 4(a) is of 0{1/N)^, because the summation with respect 
to both ko- and a is performed as seen in eq. (2.34). On the other hand, those of Fig. 4(b), Fig. 
4(c) and Fig. 4(d) are of 0{1/N)^ because one summation over remains. For example, in 
Fig. 4(b), the summation {I/Nl)^ is performed, so that this diagram gives the order 

of (l/iVL)Efc' - 0{l/N). Similarly, in Fig. 4(c), the summation (I/A^l)^ Ei E.' Efc' Efc" 
is performed. This also gives {I/Nl)"^ Eo-' Ea,-' Efc" ~ (1/^l) Efe" ~ 0{1/N). Then, for the 

a' a' a' 

leading order in 1/A^, the Dyson equation for the conduction electron, the slave-boson, and 
pseudo-fermion Green function are given by the Feynman diagrams illustrated in Fig. 5. It 
is noted that S^p vanishes at this order of approximation. Hereafter, we discuss the problem 
within the accuracy to the leading order of 1/A^, although there is no difficulty in taking into 
account higher order corrections in 1/A^. 
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Fig. 5. Diagrammatic representation of Dyson equations for single particle Green function within ac- 
curacy of (l/A^)". Thin solid, wavy, and dashed lines represent the Green function of the conduction 
electron, slave-boson, and pseudo-fermion, respectively. Thick lines denote the renormalized Green 
function for each particle. 



2.4 Dyson equation to the leading order in (l/N)^ 

When treating the problem under zero magnetic field or without magnetic order, we 
can neglect the off-diagonal self-energy part, with respect to the spin, of the conduction 
electrons, because of the properties of the Clebsch-Gordan coefficient. Indeed, from the relation 
of spherical harmonic function {Y'j'"(J7fc)}* = (— l)™YJ~'"(r2fc), we can find 

5]Kr<x^kr,j = 0, (2.35) 

r± 

J^|T4raP = J^l^krd', (2.36) 

r± r± 

where a = —a and X]r± ™eans the summation with respect to the Kramers doublets char- 
acterized with |r) = I lb 1), I lb 2), and | ib 3). For F = | ib 1), eq. (2.35) means that 
^k|+i>o-^k|+i)a- + ^k|_i)CT^k|-i)CT = is Satisfied. The analytic form of the relevant self-energy 
parts are given by 



Ska(iwn)= Ihji \-y^\Vkra\^TS2Bi{wn)F^^iiuJn + ii^n)/{Qi)>], (2.37) 
r± u„ 

ni(iz^n) = ^Y.Y.Y1 \Vkra\^TY,Gka{iuJn)F^ri^^n + iz^n), (2.38) 

{Q^)x = Y^{nf^^)x + {nbi)x■ (2.39) 
r± 

Here, {hfir)\ = {fJrfiT)x and {fn,i)x = {h\hi)x are given by the diagrams shown in Fig. 6 
whose analytic forms are given as follows: 

(n/.r)A=T5]i^°r(i^n) (2.40) 

- ^ J]5^|Vkr<x|'T2^ J]F0r(ic^„ + iM„)2Gk.(ic^„)i?i(iz^n), (2.41) 
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{nhi)x = -TY,B^{ii^n)- (2.42) 

Taking the summation with respect to the Matsubara frequency in eqs. (2.37) and (2.38), and 
performing the analytic continuation ia;„ — > a;+ = uj + iO+ , we obtain the imaginary part and 
real part of the self-energies as follows: 

ImSk.(cu+)= lim y |ykra|'fl + eMe-^(^»+^'r) 
r± 

X ImBii-uj + Xi + EiT + iO+)/(Qi)A, (2.43) 

1 1 

ReSk,(c^+) =--P — — ImEk<x(w'+), (2.44) 

^ J-oo -Uj' + lU+ 

Imni(c^+) l^kr.l'/(-'^ + + ^*r) 

^ T± a 

X ImGk<x(-^ + A, + ^ir + iO+), (2.45) 

ReU,{LO+)=--P dio' ——lmlli{uj'+), (2.46) 

71" J-oo u; - + iO+ 

where /(w) = {e^'^ + 1)^^ is the Fermi distribution function (/? = 1/T). The real part of the 

self-energies are given by the Kramers-Kronig relation from the corresponding imaginary part 

of the self-energy. Of course, we can calculate Sko- and Ilj analytically without separating the 

real and imaginary parts. However, we express them in such a way for convenience of the 

numerical calculation. 

Here, we introduce characteristic parameters £"0 and a. From eqs. (2.17), (2.45), and (2.46), 
we can find the following relation at low temperature (T < Eq/10, as shown later). 

-hmBiiio+) =a{T) ■S{uj + Xi + ef + Eo{T)) + C{uj + Xi + Sf), (2.47) 

where C{u> + Xi+ef) is a broad continuous function that is non- vanishing only in the region of 
a; > 0, and the binding energy Eq{T) and the residue a(T) of the slave-boson are determined 
through the coupled relations 

EoiT) =£f- Ren,(Ai + £/ - Eo{T)), (2.48) 

^ = 1 - ^(ReH,(u;))U=A.+e,-i?o(T)- (2.49) 

The parameter Eq = i?o(0) corresponds to the Kondo temperature Tx in the case of impurity 
Anderson model. To express a series of equations more concisely in the actual numerical 
calculation, we introduce the following notations: {Qi)\ = e~l^^^^~^^f~'^°^'^^^ (Q) , B{uj) = Bi{uJ+ 
Xi + Ef) and n(a;) = Ili{uJ + Aj + £f). Then, we can rewrite eqs. (2.37)~(2.42) as follows: 

ImSk.(^H-) = Y: \y^r.\'{l + e^-)e-'^[^o(^)+^--/l 
r± 

xlmB{-uj + Eir-ef + iO+)/{Q), (2.50) 
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Fig. 6. Fcynman diagrams for (ri/ir)A, and {nu)x within the order of (l/A^)". 



X ImGka(-w + Eir-ef + iO+), 



(2.51) 



and 



EE 



-/3[£;o(T)+£;ir-e/ 



dw /(cj) — ImGk<x(w+) 

TT 



EE 



iVkr.l 

Nl 



f{u:).{—B{u: + E,r-ef) + 



I du e-'^^^o^^^+^^r-^/^^ImGk.l-^ + iO^ 
dSfw') 



dw' 



(n,,), =e-'3(^»+^/-^«(^)) X / do; e-^(-+ 



-Eo{T))-_ 



-1 



TT 



i.j'=i.j+£'ir— £/ 
-lmB{uJ+). 



(2.52) 



(2.53) 



The validity of the approximate eq. (2.47) in the low-temperature region, or eq. (2.54) below, 
for B is shown in the next subsection, while in high-temperature region (T > Eq) eq. (2.47) 
ceases to be valid. 



1 



-lmB{LO+) =a{T) ■ 5{uj + Eo{T)) + C{u;). 



TT 



(2.54) 



To extend this formulation to a finite temperature region, we must use the Dyson equation 
self-consistently (See §4). 

Finally we calculate the average number ric of conduction electrons and nj of f-electrons 
per site: 



^EE[^E' 



Nl 



(2.55) 
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With the use of eqs. (2.37) and (2.41), we obtain 



rif 



lim Y.\T^F|,iiu;n)/{Q^ 



+ ]^ E E [-^ E Gk.(ia.„)^Sk.(i^. 



Hence, the total number of electrons per site is given by 



kcr Uu^n I 



n, + nf= lim V [t V i^°r(iu;„)/(gi 



+ 



r± 
+ 



A^L ^ r ^ d(ia;„) 



At r = 0, this relation is simplified as 

1 



log G^li}^n 



^'=+^/=^EE[^E 



d(iWn' 



■log G'kJ(i'^n) 
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(2.56) 

(2.57) 
(2.58) 



(2.59) 



(2.60) 



(2.61) 



O" kc t^n 

Furthermore, we find ImSko-(0+) = at T = from eqs. (2.50) and (2.54), so that we obtain 
the total number n = nc + Uf as 

+ ^/ = ^ E E ^(-^k. - Sk.(0)). (2.62) 



n = Tic 

kff 

The total number is determined by the volume of the /c-space enclosed by the Fermi surface. 
In a word, the Landau-Luttinger sum- rule holds at T = 0^^) . 

3. Physical properties of Generalized Periodic Anderson Model at T = 



In this section, we study the properties of the generalized PAM (2.1)^^(2.4) in the limit 
of r — > 0, where eq. (2.54) (or eq. (2.47)) is a valid approximation. (See § 3.2) With the use 
of eq. (2.54), we manipulate eqs. (2.52) and (2.53) as follows. From eqs. (2.52) ~ (2.54) 

(n/^r)A =e-^(^»+^/-^«) 

-P{Eo+E,r-ef) 



o" kff 



o|H:ra| 

Nl 



-1 



do; — ImGka-(ti;+) 



vr 



EQ-Eir + Eff 



(J ka 
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(3.2) 



where a = a(0). At T = 0, exp(— /?£'o) vanishes for Eq > 0, so that 



g I Vicryl 

Nr 



duj — ImGka-(w+) 
vr 



(u;-So-^*r + e/)2^ 

-/3(A,+e/-So) 



(3.3) 
(3.4) 



{nbi)\ =ae 

Similarly, when we use eqs. (2.37), (2.43), and (2.44) at T = 0, the continuum part of the 
spectrum of the slave-boson Green function C{lo) in eq. (2.54) can also be neglected, because 
the continuum part is smaller than the contribution from a pole by a factor of exp(— /?£'o). As 
a result the self-energy of the conduction electron, eq. (2.37), and {Qi)x, (2.39), are given as 

(Qi)^ = a{l/a - i)e"'3(A.+e/-Bo) + ag-'^CA^+^z-So)^ 

= ^-m^+Sf-Eo) ^ (3.5) 

Q|Vkr<7p 
- Eq - Eir + £f ' 



(3.6) 



where we have used the relations eqs. (2.38) and (2.49). Here, we note that the self-energy 
Sko-(w) of the conduction electrons does not have an imaginary part at w = ^ = 0. Hence a 
quasi-particle band is formed and the system behaves as a Fermi liquid, so that the Landau- 
Luttinger sum-rule holds, as seen in the previous section. From eq. (2.7), we rewrite eq. (3.6) 
as follows: 



ReSk^(tj+) =^ 



-alK 



k|±3><Tl 



± 



E, 



+ 







Q|'^k|±2)<Tp 
Eq- Ai 



+ 



«l^k|±l>aP 



UJ-4 



Eq 



(3.7) 



where, ^_|_ means the summation over the Kramers doublets. We neglect the k-dependence of 
VkFcri and we assume the summation for Vkro- over each of the Kramers doublets specified by 
r to be constant, i.e., X^-i-lVkro-P = Z^±l^ro-P = V'^ for simplicity. However, this assumption 
does not qualitatively change the physical properties. Then, we obtain 

aV^ aV^ 



ReSk<x(t^+) 



.UJ+ 



Er 



+ 



+ 



Eq — Ao 



(3.8) 



1 



^kcr — Sko-(^i^) ' 



4 



Ai 



(3.9) 



kcr 
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where and defined by the pole and residue of the Green function of the conduction 
electron, must satisfy the following equations 

"L - Ska - ReSk.(4^) = 0, (3.10) 

and 

■ _ (4^ - Eo){ai^ -Eo- Ai)(4^ -Eo- A2) 

where = a^'^. The physical meaning of is the quasi-particle dispersion of the Fermi 
liquid, because the conduction electrons hybridize with the physical f-electron through the 
self-energy Sko-- Eq and a in eq. (3.8) are determined by eqs. (2.48) and (2.49). The total 
number n {= nc + Uf) per site is determined by eqs. (2.55) and (2.56). 

j=i r± a ^0 + ^ir - £/ - "ka 

4 



^ -^EEEEtit^^^^^. (-3) 



n = ne + n^ = ^EEE /KJ^i + (1 - «)• (3.14) 

j=l cr k^ 

Here, re/ is also derived from 1 — a because of local constraints eq. (2.9), where a is defined as a 
residue of the slave-boson Green function in eq. (2.54). At T = 0, the residue a is equivalent to 
the number of /^-states per site: i.e., nj = 1 — a. By solving a set of self-consistent equations 
(3.12)~(3.14), we can obtain physical quantities at T = 0. The solution derived from the 
above equations is consistent with that of slave-boson mean field theory at T = 0. In the 
following subsections, we show the numerical results. 

3. 1 Dispersion of renormalized band 

We assume that the dispersion of the bare conduction electrons is linear in any direction, 
so that the DOS of the conduction electron per spin is approximated as follows: 

pA^) ^ j;imGL(^+), (3.15) 

^ k^ 

. . jr. --^ i-D<u;<D), 
Pa{uj) = { (3.16) 
(otherwise) . 

We adopt the following parameters; D = 1, n = ric + Uf = 1.4, V'^ = 0.02D^, = 
—0.7D, where e^p is defined as the energy level of f-electrons measured from the center of the 
conduction electron band, and Sf = — /x. Although one might think that the parameter 
V has too small a value, the binding energy of the slave-boson Eq is found to be Eq ~ 0.01-D 
in the case of Ai = A2 = under these parameters. 
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Fig. 7. Schematic structure of the renormalized bands for the Fermi level set to zero. £/ is the 
atomic f-levcl measured from the Fermi level. Ai and A2 are CEF splittings (Ai < A2). 



We show the dispersion of renormalized bands a^, given by the solution of (3.10), for 
two cases of Ai and A2 in Fig. 7(a) and (b), where the Fermi level is set to be zero. In 
the absence of CEF splitting (Fig. 7(a)), the conduction electron mixes only with two linear 
combinations of 6- fold degenerate states of f-electrons, so that such a band consisting of two 
linear combinations is renormalized, and others remain un-renormalized. On the other hand, 
in the case of Fig. 7(b) where CEF splitting exists, the conduction band is renormalized by 
hybridizing with each f-level, so that the dispersion of the renormalized band is divided into 
four bands. The insets of Fig. 7(a) and (b) show that the Landau-Luttinger sum-rule holds 
because the Fermi wavenumber is fixed as kp = 0.77r corresponding to n = nc + nj = 1.4. 

We display the relationship between £"0 and in Fig. 8(a), where other parameters 
are fixed at their previously given values. One can see that the characteristic temperature £"0, 
which corresponds to the Kondo temperature in the impurity problem, increases with increase 
in the hybridization parameter. Since the effect of pressure is simulated by an increase in 
the hybridization parameter between the conduction electrons and the f-electrons, Eq rises 
monotonously with increasing applied pressure. Fig. 8(b) shows the V"^ dependence of the 
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Fig. 8. (a) Eq vs V^, and (b) a vs for Ai = A2 = 0.0_D (dotted line), and Ai = 0.02D A2 = 0.04D 
(solid line). 



residue a of the slave-boson Green function. This shows that (= 1 — a) approaches unity 
with decreasing V . This behavior is consistent with that of the renormalization factor q derived 
from a mean-field-type approximation^^-* and a Variational Monte Carlo study^^^ on the basis 
of the Gutzwiller ansatz: 

1 l-nf/2 

q-^ = i^. 3.17 

l-Uf 



3.2 Spectral weight of slave-boson 

We calculate the spectral weight of the slave-boson at T = 0. The spectral function pb{uj) 
of the slave-boson is given by 

r \ A ^ 0+-Imn(a;+) 

p,(u;) . --ImBiu^) = - . [,^,^_R,n(u;^)p^[o^_i,,n(u;^)]. - (3-18) 

We adopt following parameters; D = I, n = Uc + rif = 1.4, = 0.021)^, e^'^'' = —0.7D, and 
show the results in the absence of Ai, A2 and Ai = 0.02D, A2 = 0.04D in Fig. 9 and 10. The 
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Imn/D 




(O/D 



Fig. 9. Self-energy of the slave-boson n(w) for the case of Ai = A2 = 0, and Ai = 0.021) A2 OMD. 
Insets show magnified section around uj ~ 0. 



function of the self-energy of the slave-boson is illustrated in Fig. 9 for two sets of the 

CEF level scheme Ai = A2 = 0, and Ai = 0.021) and A2 = OMD. At T = in the absence of 
CEF splitting, 1011(61;) has a sharp structure at u; = due to the Fermi distribution function, 
i.e., from eq. (2.51) 

Imn(u;+) = j^Y.Yl /(-^)ImGk<.(-a; + iO+), (3.19) 

where we note that Eir = Sf in the absence of CEF splitting in eq. (2.51). When we take into 
account the CEF splitting, three-step-like structure appears, because Imn(u;) is given by 

Imn(..+) =^ Yl E [/(-'^)ImGk.(-a; + iO+) 

+ ^ f{-uj + AOlmGk.(-a; + A; + iO+)] . (3.20) 
1=1,2 

Fig. 10 shows the spectral function of the slave-boson Pb{^^), for the same sets of the 
CEF level scheme. One can see that the resonant peak is exhibited at u> = Eq and a broad 
continuum appears in the region of a; > 0. The condition, to + £f — Ren(u;) = 0, is satisfied 
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Fig. 10. Spectral function pb{uj) of the slave-boson. In panels (a) and (b). the region ~ is enlarged, 
and the real part of n(w) is also shown. 



only when u) = —Eq {Eq > 0) in the region of w < in Fig. 10, while Imn(a;) = at lj < 
in Fig. 9. For this reason, pb{uj) has a resonant peak at cj = —Eq, such as 1/{uj + Eq + 10+). 
In the region of < u; < Wcdgcj where Wedge is defined by a condition ImGko-('W > i^edgc) = 0, 
Imll has a finite value in Fig. 9, so that Pb{^^) has a broad peak in this region. This result is 
consistent with eq. (2.54). 



3.3 Spectral weight of conduction electron and f- electron 

Next, we calculate the one-particle spectral weight pc{oo) of the conduction electron and 
p^j{uj) of the f-electron. The spectral weight pc{uj) is defined by 



]^ E E ^ka(w) = ^ 5] 5^ [-^ImG^^iu,, 



^ j=i 

and the spectral weight of f-electrons is as follows: 



-- fim ^ ImG^'p(u;+)/ (ft) a 



vr A, 



(3.21) 



(3.22) 



r± 
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(a) 




<Q,>x 



(b) 




<Qi>x 




a' Y" ° <^ r"a' 



<Qi>x 



Fig. 11. Lowest order contribution to G'^- in 1/A^-cxpansion. 



where 



ImG}(cj+) 



e-'-(Tr[/,r(r)6j(r)6,4]);,dr, 

At ^OO 



(3.23) 

(3.24) 
(3.25) 

r± r± 
Here, it is noted that the annihilation operator /phys,j of the physical f-electron at the i-th 
site is given by the product of the operators of the pseudo-fermion and the slave-boson as 
/phys,i = firb]- Namely, the Green function for f-electron is obtained by the convolution 
of slave-boson and pseudo-fermion Green function. The lowest order contribution to Gj in 
the 1/A^-expansion is illustrated in Fig. 11. The diagrams (b), (c), (d) in Fig. 11, when the 
summations on a were not performed, give higher order corrections to the self-energy of 
the conduction electron, according to the classification rule of the Feynman diagrams in §2. 
However, when we calculate the spectral weight, an extra summation with respect to a must 
be performed, so that the terms are returned to the lowest contributions. The first term (a) 
in Fig. 11, which is denoted as Gj(u;+)(a), is given as follows: 

G}(iu;„)(a) = - hm V T V B,(izy„)i^o (ia;„ + iz.„)/(gi)A, (3.26) 



r± 



where Bi{ivri) and F^^{iuJn + ii^n) are given by eqs. (2.17) and (2.18), respectively. The con- 
volution of Bi and F^p, eq. (3.26), is equal to Ska/I^kro-P as can be seen in eq. (2.37) if the 
k-dependence of VkPo- is neglected as in the present approximation scheme. On this basis, with 
the use of eq. (3.6) for Sko-, we obtain 

G/(^+)(a)=E' 



r± 



Eq — Eir + £/ 



(3.27) 
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(3.28) 



r± 



where £"0 > 0, so that Eq + Eiy — £/ > 0. ImGj(Lij+)(a) h^-s a value only in the region of > 0. 
Therefore, the contribution from the first term (a) in Fig. 11 vanishes. The terms (c) and (d) 
in Fig. 11 also have a value only in the region of u; > 0, so that they do not contribute to 
the average number of f-electrons. As a result only the second term (b) in Fig. 11 has a value 
in the region of w < and contributes to the average number of f-electrons at T = 0. For 
this reason, when we discuss the spectral weight of the f-electron in the region of to' < at 
T = 0, we should only take into account term (b). The analytical form of the second term (b) 
at r = is given by (See Appendix A for details) 



1 



)(b) =]^EEEl^ra| 



de/(e)[— ImGk.(-e + iOH 
L vr 



X I —lmB{£' + iO+)l B{oj+ + e + e')- 



+ e' -EiT + Ef) 



-KQ) 



+ de/(e)[^ImGk.(e + iO 



-1, 



lmB{£' + iO+)l B{-oj + e + e' + iO^ 



2/(^)' (3-29) 



{e + e' - EiY + Ef) 

where we have used the previous notation: {Qi)x = e~^^^'^^^f~^°\Q) . Substituting this into 
eq. (3.23), we obtain 

=];7^EEEl^r.py" defie) J d.'e-/^(^o+^')^ImB(.' + iO+)/(Q) 

1 



PA 



X [— ImGka(-e + iO+) • —lmB{uj + e + e' + iO^ 
L vr vr 



-1 



{uj + e' - Eir + Ef)^ 

-1 -1 - 1 

H ImGka(e + iO+) lmB{-uj + e + e' + iO+) • 

vr vr [e + e'-Eir + efy 

With the use of (2.54) at T = 0, this is reduced to 



^/M=EE 



r± cr 



U — En — Eir + e 



a\Vra\ 



OJ 



Eo-Eir + Ef)'- 



J de f{e)p^{-e)C{u; + e-Eo) 



+ EE/d^ fie)pAe)C{-u; + e-Eo)- 



alVro 



(3.30) 

(3.31) 
(3.32) 
(3.33) 



r± cr 



{e-E^-EiT + efY' 

Here, po- is defined in eq. (3.15), and C{uj) is the broad continuous function in eq. (2.54). 
The contribution from eq. (3.31) corresponds to the coherent part of the spectral weight of 
f-electrons, and that from eqs. (3.32) and (3.33) give its incoherent part. 
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Fig. 12. Spectral weight of the conduction electron and f-electron. Insets are magnified views around 
w = 0. 

We show the numerical results of Pc{oj)-, eq. (3.21), and /Cj(<j-'), eqs. (3.31)~(3.33), in Fig. 
12, where we adopt the following parameters: D = 1^ n = ric + Uf = 1.4, V"^ = 0.02Z)^, 
= —0.7D. A resonance-like peak is observed at lu = Eq in p^j{u>), the spectrum of the 
f-electron. This peak arises from the term, eq. (3.32), in which C{uj) is a broad function of 
Pb{uj) in the region of a; > 0. Physically speaking, it reflects a excitation process from the 
binding state of the slave-boson to continuum states whose minimum excitation energy is 
given by Eq. This result can not be obtained by a simple quasi-particle picture of the Fermi 
liquid, where the renormalization factor z is estimated as z = (1 — dTj{uj) / du>\uj=Q)~^ , but is 
related to a subtle structure of incoherent states with excitation energy of the order Eq ~ Tp, 
Tp being a renormalized Fermi energy. Here, with the use of these spectral weight functions, 
we find that the following sum-rule holds within an error of a few percentage at worse. 

/CXD 
do; /(w) p,{u;) + p'fiu;) = n. (3.34) 
-oo 

3.4 Kondo temperature 

Eq defined as the binding energy of the slave-boson at T = corresponds to the Kondo 
temperature Tk defined in an impurity version of the Anderson model. In the impurity prob- 
lem, the conduction electron is not renormalized, so that we can put Eq — > Tk, oP^^ — > Sko- 
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Fig. 13. (a) A2 vs Eq in the case of Ai = 0. (b) A2 vs Eq in the both case of Ai = (dot hne) and 
Ai = 0.5D (soUd line). 



7^(0) 



and — > 1. Then, we only have to take into account eq. (3.12), leading to 

Similarly, we can calculate the Kondo temperature T^^ in the absence of CEF splitting. Then, 
we can derive the following relation between Tk and t|?^ under the condition Tk <C Ai ~ 

A2 <. \ef \ < \D\, 

'"'^ --Ai-A2-Tk. (3.36) 

(See Appendix B for derivation) This result is consistent with previous works on the impurity 
problem.^''' 

Next, we consider the lattice problem, for which we must solve coupled self-consistent 
equations (3.12)^^(3.14). By solving these equations, we can evaluate £"0, which is regarded 
as the Kondo temperature in the lattice case. We illustrate the numerical results for the 
effect of CEF splitting on Eq in Fig. 13, where we adopt the following parameters; D = 1, 
n = nc + nf = 1.4, V"^ = 0.02L>^ ef^ = -0.7D, and (a) Ai = O.OD, and (b) Ai = 0.5D. 



22/37 



J. Phys. Soc. Jpn. 



Full Paper 



In case (a), if A2 = 0.0 D, Eq should be in agreement with , which is defined as the 
Kondo temperature in the absence of CEF sphtting in the lattice case. Here, with the use 
of the relations (3.35) and (3.36) derived in the impurity case, we can roughly estimate E^^ 
from Substituting the present set of parameters: po = 0.5/D, Sf = e^p — fJ- = —0.23-D, 
and y2 = 0.02L>^ we obtain 



Numerical results in Fig. 13(a) show Eq ~ lO^^D, which is reasonable compared with the 
rough estimate of eq. (3.37). 

One can see in Fig. 13(a) that £^0 decreases rapidly with increasing A2, and Eq saturates 
at A2 — 4Z). This shows that the Kondo temperature shifts from that with 6-fold degeneracy 
to that with 4-fold degeneracy. Such a behavior was not derived from the SU(A^)-PAM-|-CEF 
model, in which Eq becomes a negative value at large A ~ 4D. Namely, the application of the 
1/A^-expansion method to SU(A^)-PAM fails when we take into account the CEF splitting. 

In case (b), Eq also decreases rapidly with increasing A2 as shown in Fig. 13(b). This 
also shows that the Kondo temperature shifts from that with 4-fold degeneracy to that with 
2-fold degeneracy. We find that the value of Eq ~ 4.7 x lO^^D at Ai = O.OD, A2 = 0.5D is 
consistent with that of Eq at Ai = 0.5D, A2 = O.OD. 

4. Physical properties at Finite Temperatures 

In this section, we discuss the resistivity over the entire temperature range. To this end, 
we perform the calculations beyond the low-temperature approximation given in §3. In other 
words, when we discuss the physical properties at high temperatures (T ^ Eq), it is not 
appropriate to use eq. (2.54) given in §3, since Imll has the width of ~ T near u; = due 
to the thermal smearing effect of the Fermi distribution function. Namely, the picture of the 
resonant peak at Eq ceases to be valid so that we must take into account the entire structure 
of the spectrum of the slave-boson and the contribution from the pole of pseudo-fermion at 
high temperatures. Then, we must solve the Dyson equations eqs. (2.16)^^(2.18) for the single- 
particle Green functions self-consistently for all values of frequency uj at each temperature. 
Details of the self-consistent calculation are discussed in Appendix C. 

4-1 Spectral weight of conduction electron 

The spectral weight of Pc(t^) of the conduction electron, given by eq. (3.21), is shown in 
Fig. 14 for a series of temperatures. The parameters adopted are the same as for Fig. 14~Fig. 



19 : D = I, n = nc + nf = 1.4, = 0.021)^ ey' = -0.7 D. Fig. 14(a) is for the case without 
CEF splitting (Ai = A2 = 0), and Fig. 14(b) is for the case with CEF splitting Ai = 0.02L> 
and A2 = 0.041). Here, we have also assumed the DOS of bare conduction electron per spin 
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Fig. 14. Spectral "weight of conduction electron at each temperature for cases (a) Ai = A2 = O.OD 
and (b) Ai = 0.021? A2 = 0.04D. Insets arc magnified vie"ws around a; = 0. 



as follows: 

Pa(^) = TTZ^Pka = < ^ ^ ^ (4.1) 

^ [0 (otner"wise) . 

With increasing temperature, the hybridization gap in the spectral "weight of the conduction 
electron is gradually buried through pseudo-gap-like behavior. Such behavior has already been 
pointed out in a previous "work detailed in ref. 10. It is noted that, in the case of Ai = 0.02D 
and A2 = 0.041), three hybridization gaps, corresponding to each CEF level, appear in the 
lo"w-temperature limit. 

4-2 Spectral weight of slave boson 

In Fig. 15, we sho"w the spectral function of the slave-boson in the case without CEF 
splitting, i.e., Ai = A2 = O.OD. In the finite temperature calculation, Eq(T) is determined 
by the peak position of the spectral function of the slave-boson. In Fig. 15(a), we find that 
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Fig. 15. (a) Spectral weight of slave-boson pb{^)- Inset shows magnified view around lj = Q. (b) /Oh(ti^) 
for a series of temperatures. 




Fig. 16. Self-energy of slave-boson, Rcn(a') and lmY{{uj). for a series of temperatures. Inset is mag- 
nified view around lo = Q. 

£^0(4 X IQ-^D) is almost the same as the result Eq{Q) derived from the T = calculation 
in §3, -Eo(O) — 0.01. Namely, we can confirm that our calculation is performed from a low 
temperature (T < Tq) to a high temperature (T > Eq) continuously. Furthermore, we find 
that at a high temperature such as T = O.ID, the spectral weight of the slave-boson cannot 
be divided into the resonant peak and broad function like eq. (2.54). 

In Fig. 15(b), we show the results for a series of temperatures 4 x 10~^ < T/ D < 0.1. We 
find that the validity of eq. (2.54) collapses at T ~ 0.01 i?o(0). Indeed, the low-temperature 
approximation fails at temperatures higher than ii^o(0)/10. In the case with CEF-splitting, 
such a behavior does not change essentially. However, £'o(0) becomes smaller compared with 
the case of Ai = A2 = O.OL* (See Fig. 8(a)), and a three-step-like structure appears due to 
the CEF-splitting (See Fig. 9 and 10). 
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Fig. 17. Diagrammatic representation of the current-current response function. 

Next, we show the result for the self-energy of the slave-boson in Fig. 16. Imll(a;) has a 
sharp structure around ti; = 0atT = 4x W'^D, while its width broadens with increasing 
temperature due to the smearing effect of the Fermi distribution function. Therefore, Imn(u;) 
at high temperatures has a finite value even in the region of w < 0. Hence, the peak at to = Eq 
has a finite width, so that pb{uj) cannot be divided into two parts, a resonant peak (in the 
region of a; < 0) and a broad background (in the region of > 0), at high temperature 
T ~ £^0- 

4-3 Electrical resistivity 

In this subsection, the temperature dependence of the resistivity is discussed. We calculate 
the conductivity by means of the Kubo formula. Since the dispersion of f-electrons as well 
as the A;-dependence of V is neglected, the conductivity is expressed by the Kubo formula 

where fk = VkEk; K is a two-particle Green function defined as 

/•oo 

(a;+) = i / dt e-+*([4^(t)ck+q.(t),cL,(0)ck'_q.'(0)]). (4.3) 
Jo 

Since we are interested in the qualitative behavior of the resistivity, we neglect the vertex cor- 
rection corresponding to the inverse collision process in the Landau-Boltzmann equation. How- 
ever, the correction of the current vertex is taken into account through the Ward-Pitaevskii 
identity in the k-limit although this gives no net effect to the current vertex provided that 
the wave- number dependence of the self-energy Ej(k, tj) of the f-electron can be neglected 
compared with its frequency dependence. Indeed, the Ward-Pitaevskii identity for the current 
vertex^^) is given by 

i/r!:,,„,(p,,)|{o^(<,)).(0 = -;^, (4.4) 

where a' is the quasi-particle weight, and G is a one-particle Green function derived from the 
relation between the two-particle Green function and the vertex part. The right-hand side of 
eq. (4.4) can be approximated by —p/m if the k-dependence of Ej is neglected. By assuming 
implicitly that the Umklapp process works to violate the conservation of lattice momentum, 
the effect of collision is taken into account only through the self energy of Gko-- Then, luiK is 



- + ^ / r 
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T/D 



T/D 



Fig. 18. (a) Temperature dependence of electrical resistivity p{T) and z^^. Inset is temperature 
dependence of Eq. (b) Temperature dependences of number of conduction electrons ric, f-electron 
n f , and n = Uc + Uf , and e / measured from the chemical potential at a finite temperature. 



given as (See Fig. 17) 

lmKg{uj+) = - 6kk> j de f{£) — ilmGka(e + iO+) 

X Im|Gk+qa(e + + iO+) - Gk+qa(e - w + iO+)|. 

Substituting this into eq. (4.2) and performing some rearrangements, we obtain 

5/(e)xr 12 



(J 



ttNl 



EE 



de 



de 



ImGka(e + iO+) 



By omitting some constant factors, we define the reduced conductivity. 



a 



de 



dm 

de 



ImGka (e + iO^ 



(4.5) 



(4.6) 



(4.7) 



With the use of the renormalized Green function for ImGko- derived from §4.1, a is calculated 
numerically, and the reduced resistivity is defined as p = a~^. A result of the resistivity p 
is displayed in Fig. 18 and 19. In the absence of CEF splitting, the resistivity has a broad 
peak near T ~ 2Eq{0) as shown in Fig. 18(a). In the high-temperature region, T > 2Eq{0), 
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the resistivity shows — logT behavior, typical of the Kondo effect. Here, we define the mass 
enhancement factor follows: 



„i _ dEk(x(cj) 
" da; 



1^=0+ 



(4.8) 



which is relate to the quasi-particle weig ht zp = (1 - dT:{^{uj)/duj 



as follows: 



1.^=0+/ 

' = l + El^kr<.|'G}p(0)2zp"\ (4.9) 
r± 

where Gjp(u;-|_) is defined in eq. (3.25). With decreasing temperature, Zc~^ begins to increase 
at around T = Eq and saturates at a temperature lower than Tq ~ £"0/10. Namely, the quasi- 
particle is formed at a temperature lower than Tq/IO. The resistivity p exhibits a gradual 
decrease in this temperature region. In the present calculation of the leading order in in 
which the inter-site correlation of f-electrons cannot be taken into account, the Fermi liquid 
behavior p oc T^(at T <^ £"0/10) is not reproduced. In order to reproduce the T^-behavior, 
we need to take into account higher order terms of 0{1/Ny as shown in ref. 14. However, the 
higher order correction is known not to appreciably affect the temperature dependence of p 
at T > Tq.^^^ In Fig. 18(b), the T-dependences of Uf, Uc, and n = nc + rif, and £f are also 
displayed. The numerical calculation is performed under the condition that n is fixed. 

In Fig. 19, the T-dependence of the resistivity p{T) in both the cases of with and without 
CEF splitting are shown. Eo{0) with CEF splitting becomes smaller than that in the absence 
of CEF splitting. In the inset of Fig. 18(a) and 19, the T-dependence of Eq{T) is shown. We 
find that Eq{A x 10~^D) ~ 0.01-D in the case of Ai = A2 = 0. This value is consistent with 
£0(0) derived from the calculation at T = in §3. 

5. Pressure Effect on Resistivity 

In this section, we discuss the effect of pressure on the double-peak structure in the T- 
dependence of p{T). In heavy fermion compounds such as CeCu2Si2,^^ the double peaks merge 
into a single peak with increasing pressure. We discuss this problem on the basis of the present 
model (infinite-?/ J=5/2 generalized PAM under a tetragonal symmetry) in which the effect 
of pressure is parameterized as that of the hybridization V between conduction electrons and 
f-electrons. We adopt the following parameters; D = 1, n = -|- = 1.4, e^j^ = —0.7D, and 
we set the CEF parameters as Ai = 0.02D and A2 = 0.04D putting the case of CeCu2Si2^^^ 
in mind. We illustrate the results of the calculation for a series of hybridization parameters 
V in Fig. 20. For = g ji x IQ-'^D'^ ~ 1.0 x IQ-'^D'^, the resistivity exhibits a double- 
peak structure, for V'^ > 1.28 x the double-peak structure fades away, and a single- 
peak structure is obtained. Namely, we find the double peaks of the resistivity merge into a 
single peak when V is increased gradually, which is consistent with experiments under various 
pressures. 
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Fig. 19. (a) Temperature dependence of electrical resistivity p{T) in case of Ai = A2 = 0, and (b) 
that for series of CEF-splitting schemes. 



This result is summarized as follows. The Kondo temperature, which is related to the 
peak at lower temperatures, increases with increasing V (See Fig. 8(a)), so that one peak 
at a lower temperature shifts to a higher temperature monotonously with increasing V. On 
the other hand, in Fig. 20(b), another peak at a higher temperature shows only a slight shift 
to a lower temperature with increasing V. After that, the merged peak, which formed at 
V"^ > 1.28 X lO^^D^, shifts further to much higher temperatures with increasing V. This 
behavior describes the experimental results fairly well. For a hybridization V'^ = 7.2 x W^^D^ 
smaller than V'^ = 8.1 x W~^D^, the double-peak structure of the resistivity also becomes 
invisible. The parameter space in which p(T) exhibits the double-peak structure is limited, 
and is dependent on the energy scale of Ai and A2 etc. 
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Fig. 20. (a) Temperature dependence of resistivity for a series of hybridization. Arrows indicate the 
temperatures T™<^'^ and Tj^^x .^^j^^j-g _dp/dT is zero, (b) vs V in the case of the CEF spUtting 
Ai = 0.02D, As = OMD. 



6. Conclusions and Discussions 

We have studied the behavior of the electrical resistivity p{T) under pressure of Ce- 
based heavy fermion systems in the framework of a generalized periodic Anderson model 
(generalized PAM) in a manifold of J = 5/2 under a tetragonal configuration. We applied the 
1/A^-expansion method to the generalized PAM, and presented numerical results within an 
accuracy of (l/N)^ for the calculation at finite temperatures as well as T = 0. Furthermore, 
by using the present extension of the (l/A^)-expansion method, we investigated the effect of 
pressure on the resistivity p{T) by parameterizing the effect of the pressure as a variation of 
the hybridization between the conduction electrons and f-electrons. As a result, we reproduced 
the behavior that the double-peak structure in p{T) is shown to merge into a single peak with 
increasing pressure, which is in good agreement with experiments for numerous heavy fermion 
compounds. 

Here, we point out some remaining problems and perspectives. The present model may be 
too simplified for quantitative study in the sense that we have neglected the k-dependence of 
the self-energy as well as the hybridization VkPo-- We may be able to discuss the k-dependence 
of VkTa qualitatively by parameterizing the weight of the hybridization parameter V for each 
level r without setting V constant. On the other hand, the k-dependence of the self-energy 
is treated by taking into account inter-site correlation which is observed in the higher order 
terms of 0(1/A^). By studying the problem up to 0(1/A^)^, we can discuss the inter-site 
correlation such as the RKKY interaction as detailed in ref. 25. Fundamental problems such 
as the competition between the RKKY interaction and the Kondo effect around quantum 
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Appendix A: Calculation of G^(ct;)(b) 

In this appendix, we discuss the structure of Gj(w)(b) in detail. At T = 0, the analytical 
expression of eq. (3.29) is derived, while at a finite temperature T > Eq an additional term 
becomes important. From the diagrammatic representation in Fig. 11(b), we write Gj(ia;„)(b) 
as follows: 

G}(ia;„)(b) ~^^^T^TY,\yTa\'GU-OB^{il^n) 



n2 



X Fjp(ia;„ + if„)5i(iw„ + ivn + iw^)/ (Q) 



A- 



(A.l) 



Here, we take the summation with respect to the Matsubara frequency. 



G'/(iw„)(b) 



1 



27ri 



r± cr k,. 



E E E l^r^l' / [^ImGk.(-e + iO+ 



de'n{e')Bi{e')Fli\ujn + e')Bi{iuJn + e + e')/{Q. 



i X, 



(A.2) 



c 



where n(e) is the Bose distribution function. Here, we also use the notation: B{uj) = Bi{io + 
Xi + £f) and {Qi)\ = e~^^^*~^'^f~'^°'^ {Q). Then, the above equation can be rewritten as 

— ImGka(-e + i04 
vr 



C/Mcb) =- ^EEEl^r.l^ / de f{e) 



r± (T k„ 



- / de'e-^(^«+^')5(e')i^'r(i^n + ^')^(i^n + e + e')/(Q), 

TTl Jc 



27ri 
1 



r± (T k„ 



E E E l^ral' / de /(e) [— ImGk.(-e + iO+) 



/ de' e-('^'''>+'">B{e')Fl{iu;n + e')B{iu;n + e + e')/{Q), 



Nl 



r± <T k, 



E E E l^r.l' / de /(e) [--ImGk<.(-e + iO+)] /{Q) 



— ImS(e' + iO+)j Biiujn + £ + 



+ / de' e-'^(^o""+"') — Im5(e' + iO+) B{-iLOn - e + e') — 
7 L vr J e' 



-1 



+ 



de' 



e-/3(So+^')^(e')^(£ + e' + icu„)' 



:'=-iuJn+Eir-ef 



{£' -e-Eir+EfY 

(A.3) 
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Fig. 21. p^(a;)(b) vs w for a series of temperatures. 



After performing the analytic continuation iu;^ — > = u;+iO+, and changing the variable e 
— e in the second term, we obtain eq. (3.29). At T = 0, the third term, which is the contribution 
from the double pole of the pseudo-fermion, vanishes because of a factor exp(— /JE'o). 

At a finite temperature (T > Eq), however, the third term in eq. (A. 3) is non- vanishing. 
The contribution from in Fig- ll(t)), consists of both the pole of the slave- 

boson, and the double poles of the pseudo-fermion, Gj^'^^^(u;)(b)- Here, we 

define /3}(a;)(b) = -l/7rImG}(u;)(b), p}^^^^(u;)(b) = -l/7rImG}^''^V)(b), and /9}^*^(u;)(b) = 



-l/7rImG}^'^^^(a;)(b). Then, 



p}(^)(b) = p}^'''^H(b)+p}^*^M(b), 



where 



Pf'\-^)i^) = EEl^r.P / d. m [ ds'e-/^(^«+^')^ImB(.' + iO+)/(Q) 



r± o- 



-e) Im5(u; + e + e' + iO+) • 



+ p,(e) ■ —ImBi-iv + e + e' + iO+) • 

TT (e + e' - Ajr + £/) 



and 



p}('^^(u;)(b) =EEl^r.py" de /(eK(-e) • (-1; 



• e 



-/3(£;o+S,r-£/) 



7(0) 



r± o- 



[-^Im^(-w + Eir -ef + iO+)B{e + Eir -Sf) 



+ B{e + Eir - ef )lm 



dB{e' + iO^ 



de' 



e'=~uj+Eir~ef 



(A.4) 



(A.5) 
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In Fig. 21, we show the temperature dependence of p^j{uj){b) as calculated numerically. 

Appendix B: Kondo temperature in impurity case 

In this appendix, we derive the relation between the Kondo temperature and the CEF- 
splitting in the impurity Anderson model, eq. (3.36). Eq is defined as the binding energy 
of the slave-boson. At T = 0, Eq corresponds to the Kondo temperature Tk defined in an 
impurity version of the Anderson model. In the impurity problem, the conduction electron is 
not renormalized, so we can put £^0 Tk, oP^^ — > eko- and yl^ 1, in eq. (3.12): 
rj. ^ ^ /(ek<x)|VraP 



^L^^^TK+E.r-Sf- Sko 

Nr 2^ ^ Tr^ - Nt 2-^ 2-^ 



Nl^^Tk- Eka 
ko- 



Nl^^ Tk + Ai - 

<^ ko- 



y2 

Nr 



, . . , (B.l) 

Nl ^ Tk + A2 - Ska ^ ^ 

O- kcr 

where we have used the following approximation: Yl± l^ro-P = l^rsP = V^, where ^_|_ 
means the summation over the Kramers doublets. Here, we assume Tk <^ Ai A2 <^ |e/| < 
\D\, which is a reasonable assumption in an actual system, and i^/NL)Ylk is replaced by the 
integral po de, po being the DOS of the conduction electron per spin at Fermi level. Then, 
eq. (B.l) is reduced to 

+ 2V'po r de f{e)—^. (B.2) 

J-D A2 - £ 



2V^Po 



D 



■ e 
1 



de f{e) 
D i^i- e 



On the basis of this equation at T = 0, the Kondo temperature Tk is evaluated as follows: 
e-integration in (B.2) is performed explicitly, leading to 



£/ =2V polog 



^-2V'pQ log 



D-Tk 



Tk 



+ 2yVolog 



D - Ai 



Ai 



+ 2V^pQ log 



L» - A, 



D 



Tk 





D 




D 


+ 2V^po log 


Ai 


+ 2V^po log 


A2 



Then, it is easy to find a solution for Tk as 



(B.3) 



(B.4) 



Similarly, we can calculate the Kondo temperature t'^ in the absence of CEF splitting. By 
putting Ai = A2 = in eq. (B.l), we obtain 

V"^^ /(ek.) 



Nl 



EE 

cr k„ 



T 



(0) 



(B.5) 



K 



eko- 
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Hence, 



From the relations between Tk and Tj^\ we can derive equation (3.36): 



K 

3 



Ai • A2 • Tx. (B.7) 



In the leading order, the result with use of 1/A^-expansion method in case of the impurity 
problem, is consistent with the previous impurity study. ^'''^^^ 

Appendix C: The algorithm for self-consistent l/AT-expansion 

In this appendix, we explain the algorithm of self-consistent 1/A^-expansion over the entire 
temperature range. The whole procedure is summarized as follows: 
1) Self-energy of conduction electron 

First, we choose the self-energy Sko-C^i^) of the conduction electron as a trial function. If we 
start from the low-temperature region, we can choose Skcr(w) derived from low-temperature 
approximation, eq. (3.6), i.e., 

where Eq and a is obtained by the coupled self-consistent equations at T = 0, eqs. (3.12) 
~ (3.14). We restrict the range of the energy, from — Wmax to Wmax- In an actual numerical 
calculation, lo is defined as a discrete value ujj {j = -jma.^, -Jmax + 1, • • • , jmax - 1, jmax), 
where jmax is an integer of maximum cut number of the energy. For example, we define ujj as 
follows: 

= Wmax- (C.2) 



Jmax 

We use Wmax = 2D and jmax = 700 in the present paper. 
2) Green function of the conduction electron 

The Green function of the conduction electron is obtained from the above self-energy as 

G^..{^) = ^ , y (C.3) 

where ^ is the chemical potential. Then, {^) Yl\i„ G'kcr(i^) is performed as follows, 

1 /"^ 
^ VG'ka(w) =Po / de- 



cj-e + /i-Ska(w)' 

^ + + ^ka(^) ,^ .s 

2D ^uj-D + fi- Sk<,(u;) ' ^ ' 

where po = 1/2D as in eq. (3.16), and /i = in our numerical calculation. 
3) Self-energy of slave-boson 

The imaginary part of the self-energy ImII(u;), eq. (2.51), of the slave-boson is calculated 
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[Initial Parameters] 

total number : n (fixed) , atomic f-level : Ej (fixed) , hybridization parameter : V ' (fixed) 
binding energy of the slave-boson : Eq residue of the slave-boson : a 
chemical potential : p. , initial temperature : T 



|i is fixed 



[self-consistent IW-expansion routine] 

At a certain temperature T-^ 

DEka 2) ^ 3) n ^ 4) B ^ 5) Eq 6)<Q> 7) E^a 

t I 



n„ew= "c+nf 
\ 

|X is redefined 
next step (J; Ti^i ) . ^ 



physical quantities 



Fig. 22. An example algorithm for self-consistent l/7V-expansion. 



by using Gk^, eq. (C.3), as 

Nl 



Imn(a;) =Y,Y.Y. ^^/(-^ + - ef)lmG^,{-u + Eir - Sf). (C.5) 



r± <T k,, 

Ren(a;) is calculated by the Kramers-Kronig relation. Hence, we can obtain n(u;) = Ren(u;)-|- 
i • Imn(a;) 

4) Green function of slave-boson 

The Green function B{(jj) of the slave-boson is given as 

uj + £f — llyoj) 

5) Determination of Eq{T) 

The binding energy of the slave-boson Eq{T) defined as the real part of the pole of B{uj), 
eq. (C.6). However, in actual calculation, we approximate it as the value which satisfies the 
following equation: 

Sf - Eo{T) = ReUi-EoiT)). (C.7) 

This equation is an approximation of eq. (2.54). 
6) Calculation of {Qi)\ 

We have defined {Qi)x = e~^^^^~^^f~^°^'^^^ (Q) , where Eq{T) is obtained by procedure 5). 
We calculate {Q) from eqs. (2.39), (2.48), and (2.49) by using the Green function of the 
conduction electron in procedure 2), the slave-boson in procedure 4), and the pseudo-fermion. 
7) Self-energy of conduction electron 
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Finally, through the procedures 2) ~ 6), we can calculate Sko- again by eq. (2.43) as 

Im^U^) = |Vkr<x|'(l + e^-)e-^(^o(^)+^--e,) . + E,r - Sf). (C.8) 

r± 

ReEko-(w) is derived from the Kramers-Kronig relation. As a result, we can obtain the newly 
defined Skcr('^) = ReSko-l'^) + i • ImSko-l"^)- 

This series of procedures, 1)~7), forms a self-consistent cycle. We show the concept chart 
of this cycle in Fig. 22. In our calculation, this procedure was performed with fixed n, the 
total number of electron per site. In order to check the self-consistency, we define the following 
amount, 

Jmax 

|Imnnew(wj) - Imnoid(wj)| 
Sluiti = tzi^ . . (C.9) 

2jmax ~l~ I 

ImIIncw('^) is a new Imll((x;) obtained from one cycle by using Imnoid(w). This self-consistent 
cycle is repeated until 5Imn(a;) becomes less than 1Q~^D ~ 10~^-D. By performing this 
procedure for each temperature, we can calculate the temperature dependence of physical 
quantities such as the electrical resistivity. 
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